Separation-dependent localization in a two-impurity spin-boson model 
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Using a variational approach we investigate the delocalized to localized crossover in the ground 
state of an Ohmic two-impurity spin-boson model, describing two otherwise non-interacting spins 
coupled to a common bosonic environment. We show that a competition between an environment- 
induced Ising spin interaction and externally applied fields leads to variations in the system-bath 
coupling strength, ac, at which the delocalized-localized crossover occurs. Specifically, the crossover 
regime lies between Qc = 0.5 and = 1 depending upon the spin separation and the strength 
of the transverse tunneling field. This is in contrast to the analogous single spin case, for which 
the crossover occurs (in the scaling limit) at fixed etc f« 1. We also discuss links between the 
two-impurity spin-boson model and a dissipative two-spin transverse Ising model, showing that 
the latter possesses the same qualitative features as the Ising strength is varied. Finally, we show 
that signatures of the crossover may be observed in single impurity observables, as well as in the 
behaviour of the system-environment entanglement. 



I. INTRODUCTION 

Any real quantum system is surrounded by some un- 
controllable environment, interactions with which gen- 
erally lead to differing behaviour from that expected if 
the system were somehow isolated The spin-boson 
model [5, '3'] is a popular starting point for investigations 
into such dynamics, encapsulating the effects of quantum 
decoherence, dissipation and relaxation on the otherwise 
coherent spin evolution. Furthermore^he model exhibits 
non-trivial ground-state behaviour [2|-l8|, displaying a 
zero-temperature (quantum) phase transition [3, as a 
function of system-bath coupling strength, attributed to 
zero-point (rather than thermal) fluctuations within the 
bath. Besides being of general theoretical interest, many 
physical systems in the solid state, and elsewhere, are well 
described by models of a spin-boson type 0, H, [Tll - fTsj . 
Specific experimentally relevant examples include large 
arrays of trapped ions [14l] . the persistent current in a 
metal ring threaded by an Aharonov-Bohm flux [l^ [l^ , 
and atomic dots coupled to a Bose-Einstein condensate 
bath ^17,] . These systems are of particular importance, 
since it is predicted that they show qualitative and de- 
tectable changes in ground-state properties as a function 
of accessible external parameters. 

The spin-boson model considers a two-level quantum 
system, such as a spin-1/2 particle or a magnetic impu- 
rity, interacting with an (infinite) bath of harmonic oscil- 
lators representing the environment. The corresponding 
Hamiltonian is generally written in the form (for h ~ I) 

H = ;^CT^-^cr^ + o-^ V(gk&k + ffk^k) + Vt^k^lL^k, (1) 
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where e is the energy bias between the system states, 
A is the (bare) tunneling strength, and <Ti {i = x, y, z) 
is the usual ith-Pauli operator in a basis where cr^ = 
|0)(0| — The bath is represented by the creation 

(anihilation) operators b^^ (6k) for each bath mode, with 
wave- vector k and corresponding angular frequency Wk- 
The system-bath interaction is captured by the coupling 
constants g^. 

As is well known, the interaction of a quantum system 
with an environment of the type given in Eq. ([Ij causes a 
renormalization of the bare system energy levels and, in 
particular, a suppression of any tunneling probability the 
system may possess [ll-Q . For the spin-boson model, the 
system-bath interaction can be completely characterized 
by the spectral density J{uj) = J2k l5kP<^('^ ~ ^^k), which 
we shall take here to be of the paradigmatic Ohmic form 
J{id) — {a/2)uj for uj < lJc, where a is a dimensionless 
coupling strength and ujc a high- frequency cutoff [H-Q ■ In 
this case, it has been found that above a certain critical 
system-bath coupling strength, ac, the tunjneling proba- 
bility is completely suppressed (A — )• 0) For small 
e/A, as the parameter a is increased through the 
ground state of the two-level system shows a crossover 
from being dominated by the tunneling term {A/2)ax, 
and hence delocalised, to being dominated by the bias 
term (e/2)a^,, and therefore localized in either |0) or 
|1) [1, [iMl- At zero temperature and for e = this local- 
ization phenomenon has been identified as a Kosterlitz- 
Thouless quantum (rather than thermodynamic) phase 
transition [l^ . Calculation of the Ohmic critical cou- 
pling strength in this regime has found ac ~ I for small 
A/ujc i. 

In this paper, we shall investigate the delocalized- 
localized crossover in the ground state of a pair of non- 
interacting two-level systems in a common bath of har- 
monic oscillators, termed here the two- impurity spin- 
boson model. In particular, we elucidate how this 
crossover depends upon the separation between the impu- 
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rities through a bath-induced inter-spin interaction. Aside 
from being a natural extension of the single impurity 
model, the two-impurity case represents perhaps the sim- 
plest dissipative model in which to explore the interplay 
of coherent system interactions and the dissipative in- 
fluence of the bath. This has relevance, for example, in 
the field of quantum computation [20|, where the two- 
impurity model could be thought of as the basic unit of 
a dissipative spin chain l2]]42a or as two quantum bits 
in a dissipative register [ill . 124 [25| . The model has also 
recently gained attention since it allows for the study of 
bath generated correlations and entanglement shared be- 
tween the impurities (26l - [30{ . 

To perform our analysis we employ a variational tech- 
nique originally developed by Silbey and Harris [3]| . 
The method consists of assuming a particular variational 
form of the ground-state wavefunction of the combined 
system and bath, and a subsequent optimization based 
upon a minimization of the associated free energy (or 
ground state energy at zero temperature). While this 
technique might be vulnerable to errors in certain lim- 
its, it has proven to be relatively robust when applied to 
single spin-boson systems described by Ohmic spectral 
densities [3]| . Furthermore, conclusions drawn from the 
method have also been verified by path integral 0, [s^l , 
flow equation [s^ and scaling techniques [34 |. as well 
as by Bethe-ansatz [l^, [s^, [Sal, Numerical Renormal- 
ization Group [l^. Density Matrix Renormalization 
Group 0, exact diagonalization [s^], and Monte Carlo 
calculations [s^ . Besides its relative simplicity, the varia- 
tional technique is also attractive from the point of view 
of gaining insight into the form of the ground-state of 
the model, and how this varies through the delocalized 
to localized crossover. Furthermore, as we shall show be- 
low, it can be used to provide analytical calculations of 
bath-induced spin interaction terms and tunneling renor- 
malizations inherent to the model. 

While it is known that a single un-biased two-level sys- 
tem immersed in a bosonic bath will enter its localized 
regime as the system-bath coupling strength a — >• 1 (for 
A/u!c ^ 1), it is not clear whether the same conclusion 
holds true for a pair of two-level systems immersed in 
a common bath. In this case, it is possible that the 
delocalized-localized crossover could occur at a differ- 
ent point due to an induced, bath-mediated interaction 
between the spins that alters the system energy struc- 
ture [26l - [29| . For the model considered here, the mutual 
interaction with the bath induces an Ising-like coupling 
between the impurities with a strength that depends 
upon their separation. We find that for closely spaced 
impurities, corresponding to a strong Ising strength, the 
crossover region occurs around ac = 0.5. As the spin sep- 
aration is increased, the Ising strength decreases, and we 
find Qfc — ^ 1 as the impurity separation becomes infinite. 

The paper is organized as follows. In Section [ll] we 
review the application of the variational approach to the 
ground state of the standard single impurity spin-boson 
problem fsj, demonstrating its validity in comparison 



to a Bethe-ansatz calculation [1,[3|. In Section Uni we 
define the two-impurity model and apply the variational 
transformation in this situation, elucidating the origin of 
the bath- induced Ising interaction. In Section HVl we em- 
ploy an approximation on the induced Ising term that 
allows a straightforward identification of the delocalized- 
localized crossover region in a number of regimes, while 
in Section |V] we perform the full variational calculation 
without such a simplification. In Section lVTl we study the 
variational ground state in more detail and propose sig- 
natures of the crossover behaviour, while in Section IVIII 
we summarise our results. 



II. SINGLE IMPURITY SPIN-BOSON MODEL 

Before we go on to study a pair of two-level systems in 
a common bath, it is instructive to apply the variational 
technique to the (single) spin-boson model, Eq. ([1]), in 
an effort both to understand the variational method em- 
ployed, and also to assess its validity. 

Let us start by considering the ground state of the 
Hamiltonian with e = 0: 



:J2^9kbl + g*M + J2^i^blbk. (2) 



In the limit — )• with A ^ 0, the spin is entirely de- 
coupled from the bath and its ground state will be the <Tx 
eigenstate (l/\/2)(|0) + |1)). The state of the bath will 
be some superposition of number states (eigenstates of 
6jj.6k) that depends upon the temperature. In the oppo- 
site limit, A -> with ^ 0, the system-bath interac- 
tion now dominates and the oscillators constituting the 
bath will be displaced from their equilibrium positions to 
minimize the corresponding interaction energy. We may 
write the ground state of the combined system-plus-bath 
in this case as 



!*> = 



^o) |o) 



(3) 

where \Bo) is the ground state of the bath for vanishing 
system-bath coupling, and we have defined the displace- 
ment operators [39] 



exp 



± 



/5k 



(4) 



In the general case, when neither limit is met, the 
spin-boson Hamiltonian is not straightforwardly diago- 
nalizable. Note, however, that in the state described by 
Eq. ^ each oscillator is displaced by an amount deter- 
mined by the the ratio gk/wk, and that as grk — >■ the 
correct uncoupled ground state is recovered. The varia- 
tional theory thus assumes that the ground state of the 
spin-boson Hamiltonian for non-zero gk and A is always 
of the form of Eq. (jS]), but allows for the possibility that 
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the amount a given mode is displaced may have a more 
comphcated dependence on the Hamiltonian parameters. 

With these considerations in mind, we now reintroduce 
the energy bias between the spin states and proceed by 
writing down the total Hamiltonian, Eq. ([T]), in a basis 
|0) , \B+) |1)}, where |B±) = nk^(±/k/^k) |i?o), 
and we have introduced the as yet to be determined varia- 
tional parameters /k. These will be found by minimising 
the free energy of the total system [3l|. At zero temper- 
ature, we obtain 



e Ar 



(5) 



where in the new basis dz 
\B,)\l)(l\(B\ and 



|S->|0> (0|(B+| 



^=E^kVk(/k-23k). 



(6) 



Importantly, the tunneling matrix element has now been 
renormalized due to the system-bath interaction: = 
A{B), where 



(B) = (B±|S±) =exp 



2E(/kM)' 



(7) 



Diagonalization of H in the transformed basis then gives 
a ground state energy of Aq = ^(2^? — fy), where rj = 
•y/e^ + A^, and the corresponding ground state 



\(l>o) = nQ 



|i?_>|0) + |i?+>|i; 



(8) 



where no = {{il ^ £)^/^r + 1)^^^^ is a normalisation 
factor. 

The task now is to find the variational parameters /k, 
which in turn will allow us to evaluate the renormalized 
tunneling strength A^, and hence obtain the approxi- 
mate ground state. To do so, we naturally impose the 
condition that the ground-state energy, Aq, should be 
minimised. This leads straightforwardly to 



A 



(9) 



and our expression for the renormalized tunneling 
strength becomes 



Ar — A exp 



9l 



(^k + A^/r^y 



(10) 



We now take the continuum limit to convert the summa- 
tion over k into an integral with respect to uj, and recall 
the definition of the (Ohmic) system-bath spectral den- 
sity, J{uj) = X^k IffkP^Ci^ — Wk) — {a/2)uj. With these 
replacements, we find 



Ar — A exp 



{lj + Al/-qf 



(11) 




FIG. 1: (Color online) Expectation value of (Tz for the varia- 
tionally determined spin-boson ground state (plot points) and 
using Bethe-ansatz techniques (solid lines) as a function of a, 
plotted for various values of e (in units of ujc). Blue circles: 
e — 0.005; red squares: e = 10~^; yellow diamonds: e — 10"*; 
green upright triangles: e = 10~^; and purple inverted trian- 
gles: e = 10"**. In all cases A/cjc = 0.01. 



Note that had we written the original Hamiltonian in a 
basis defined with the displacement operators of Eq. ^ , 
i.e. functions of rather than /k, the integral in Eq. pT|) 
would suffer from an infra-red divergence, and we would 
conclude (incorrectly) that A^ = for all values of a. 

In the present case, the integration can be performed 
straightforwardly and leads to the following equation 
which one must solve self-consistcntly for A^: 



A, 



A2 



Al 



exp 



-auJcT] 



UJcV 



(12) 



where A^ takes on values between A and as a is in- 
creased from zero. For A/uc <C 1 and e = 0, Eq. 
gives the well-known behaviour A^ ~ A(A/a;c)"/^^~"'' [3, 

In order to assess the validity of the variational tech- 
nique, in Fig. [T] we plot the ground state magnetisation, 
(cz) = (0o| \4'o) = — e/?y, as a function of a for various 
values of e, where we have set A/uic = 0.01 and A^ has 
been found by numerically solving Eq. (|12p . Shown also 
are the corresponding plots generated by mapping the 
spin-boson model to the Kondo model and using Bethe- 
ansatz solution techniques, details of which can be found 
in Refs. [1, [H, HI]. For all values of e the methods show 
good qualitative agreement. Most importantly, the varia- 
tional calculation correctly identifies the region of a over 
which the ground state becomes dominated by the bias 
(localization) rather than the tunneling (delocalization), 
though it should be noted that (cr^) reaches its minimum 
value {{(7z) — > — 1 as Ar — >■ 0) somewhat more sharply 
than in the Bethe-ansatz calculations. We can therefore 
be confident that the variational method does capture 
the localization crossover in the ground state behaviour 
in which we are interested. 
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III. TWO-IMPURITY SPIN-BOSON MODEL 

We now return our attention to the main subject 
of this work, determining the ground state behaviour 
of a pair of impurity spins interacting with a common 
bosonic bath. Since we expect the bath to mediate a 
separation-dependent coherent interaction between the 
spins [HI, [2^, [H, m, we make their spatial separation 
exphcit by placing them at positions ri and Y2- The total 
Hamiltonian we consider is then given by 



H = 



k 



6ke 



(13) 



where tr" (n = 1, 2; i = x, ?/, z) is now the usual ith-Pauli 
operator acting on the relevant spin, and we have as- 
sumed that the system-bath coupling constants for each 
spin differ only in a position-dependent phase factor. For 
simplicity, we now limit our analysis to the case in which 
there is no bias on either spin. 

We proceed in a slightly different manner to Section HIl 
and apply a unitary transformation to the Hamiltonian 
which displaces each oscillator by an amount parameter- 
ized by the variational parameter /k. We note, however, 
that at zero temperature this procedure gives the same 
ground state as one would obtain following the method 
of the previous section. The transformed Hamiltonian is 
written H = qSi+s^ fj^-iS^+S^) ^ jj^ ^ jj^ ^^^^i 



exp 



±a;'^ak(6le'''-'"-6ke-'''-'") 

k 



(14) 



where ak — /k/'^k is assumed real. The transformation 
is aided by the observation that, provided the dispersion 
relation of the bath is isotropic and the variational pa- 
rameters satisfy /k = /-k, the commutator 



52] - 2i<jlal J2 "k sin (k • (ri - r^)) 



(15) 



where the inverse temperature is /3 = l/fc^T, while the 
form of Hs has been chosen such that the expectation 
value of Hi with respect to Ho vanishes. We shall see 
that this significantly simplifies the determination of the 
{/k} below. 

There then remain two forms of system-bath interac- 
tion. Hi = Hz + H±, where 



(18) 



and 



H, 



A 



^ ( (i?i"'-(i?»< + (i3i"'-(i?))a!! ), (19) 



with (tJ = (l/2)(cr" ± icTy ), and bath operators again 
given by products of the displacement operators: 



B 



in) 



exp 



±2^ak(6[e'''--"-6ke- 



(20) 



Note that if we assume the bath to be in thermal equi- 
librium, these four bath operators all have same ex- 



(B. 



± 'Ho 



pectation value with respect to Hq: 
exp [-2 J2k a^coth(/3wk/2)] . 

The unperturbed Hamiltonian, Hq, has two important 
features. Firstly, the tunneling strength, A^, has been 
renormalized. Secondly, the two spins are now coupled 
via a bath mediated, separation-dependent, Ising-like in- 
teraction, with a strength 



^ = E^k Vk(2.gk - /k)cos(k • (n - r2)). 



(21) 



Evaluation of both and X requires knowledge of the 
set of variational parameters {/k}- The variational pro- 
cedure determines these by free energy minimization ar- 
guments. However, before we continue the analysis, we 
outline a significant simplification which can be made. 



vanishes once the summation has been performed, re- 
gardless of the dimensionality or frequency spectrum of 
the system-bath interaction. 

The variational technique now relies on a careful choice 
of Hs, Hg and Hi from the various terms available af- 
ter the transformation. We define the new unperturbed 
Hamiltonian as Hq = Hs + Hb , with 

Hs = {al + nl) - 2Xa\Gl + 2 ^ /k(/k - 2gk)M, 

(16) 

and Hb = Hb- Here, A^. is now determined by the finite 
temperature generalization of Eq. ([7|) , and is given by 



A,. A{B) = A exp 



2E"kCoth(/3a;k/2) 

k 



, (17) 



IV. VARIATIONAL CALCULATION 
A. Crude Ising approximation 

The variational parameters {fk}, appearing in both A^ 
[Eq. (ini)] and in the induced Ising strength X [Eq. dUD], 
were introduced to overcome an infra-red divergence in 
A,, that occurs for an Ohmic spectral density when apply- 
ing a polaron transformation to our Hamiltonian, since it 
fully displaces the bath modes as in Eq. ^ [3l| . As men- 
tioned previously, this divergence would lead to a com- 
plete suppression of the tunneling probability, A^ — > 0, 
and can be seen by making the replacement /k — > .gk in 
Eq. inD, and using J{uj) = (a/2)w. 

However, we may make the replacement /k — > .gk in 
the definition of the Ising strength [Eq. ([2T|) ] and find 
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FIG. 2: (Color online) Measure of the correlation between 
the bath-induced fluctuations experienced at each impurity 
spin, plotted as a function of the scaled impurity separation 
tdOJc — |ri — r2\ijJc/c. The different curves correspond to sys- 
tem bath coupling in one dimension (blue dotted curve), two 
dimensions (dashed red curve), and three dimensions (solid 
green curve). In all cases the correlation is maximised for 
zero separation (/_d = 1, complete correlation) and tends to 
zero as the separation goes to inflnity (/_d = 0, no correla- 
tion). 



oscillations, becoming zero whenever td = mruic, for 

n = 1,2,3, .... 

. . . . . (1) 

Note that, ignoring any spatial correlations in Bj^ and 

b'^\ our transformed Hamiltonian H with the replace- 
ment X — >■ Xc now has exactly the same form as that 
which would be obtained if we transformed a Hamilto- 
nian describing two spins in separate baths, each subject 
to a transverse field of strength A, and coupled via a fer- 
romagnetic Ising field of strength 2Xc [HI 112 ■ That is, 
had we transformed the Hamiltonian 

k k 

+ ^i^k^^k&k + ^i^kflkflk, (23) 

k k 

where we have introduced a second bath which couples 
only to the second spin and is described by creation (anni- 
hilation) operators aj^ (ok), with corresponding frequen- 
cies fc'k- 



that it suffers from no such divergence. Therefore, to 
some level of approximation at least, we can make this 
replacement (in Eq. (PT|) only) and evaluate X outside 
the variational calculation. We know that in the limit 
that the coupling of the system to the bath completely 
dominates, the oscillators are fully displaced, i.e. /k — > 
gk anyway. Hence, we can identify this replacement as a 
kind of strong coupling approximation on X, as will be 
discussed in more detail in Section IVBl We shall refer to 
the Ising term evaluated within this approximation, and 
to the approximation itself, as "crude" since it does not 
take into account deviations of /k from in X. 

Assuming a linear dispersion relation |k| = w/c, where 
c is the excitation speed, we find 



Xr 



-foitdUJc), 



(22) 



where the subscript C indicates a crude value. The impu- 
rity distance dependence enters through td = |ri — r2|/c, 
which is the time bosonic excitations take to travel be- 
tween the spins, and determines the value of the func- 
tion foix) {D = 1,2,3), which is a measure of the 
(separation-dependent) correlation between the bath in- 
fluences seen at each spin, and is therefore dependent 
on the dimensionality D of the system-bath interaction. 
We find fi{x) — sinc(a;) in one dimension, f2{x) = 
1F2 ({1/2}, {1, 3/2}, -a;V4) is a generalized hyperge- 
ometic function in two dimensions, and f^{x) — S\(x)/x 
in three dimensions, Si(a;) — {sin t/t)dt being the sine 
integral function. As shown in Fig. [2l in all cases foix) 
has a maximum value of /d(0) — 1, and in two and 
three dimensions has a minimum value foioo) = 0. 
Additionally, in one dimension fi{x) displays decaying 



B. Free energy minimization 

Precisely as in the single-spin case, our task is now 
to find the set of variational parameters {/k}, which will 
then allow us to find the renormalised tunneling strength 
Ar- If, for a given a, we find that — ?> (i.e. /k — >■ 
gk), the system will be dominated by the induced Ising 
interaction, forming a ferromagnetic or antiferromagnetic 
pair. The spins will be unable tunnel between their states 
|0) and |1) and will be said to be in a localised regime. 
On the other hand, if A^ 7^ 0, the tunneling probability 
remains finite and the spins are delocalized. We expect 
that as a — ?> ac, A^ — > 0, and we enter a regime in which 
the renormalised tunneling has a negligible influence on 
the ground state. 

To find the set {/k} we follow Refs. and 31| and 
compute the Bogoliubov-Feynman upper bound on the 
free energy of the total system-plus-bath, Ab , which is 
related to the true free energy. A, via Ab > A [43], where 

Ab - ~/3-MnTr{exph/?i7o]} + {Hi)h, + ^((^Difo)- 

(24) 

We have constructed our perturbation terms Hj and sys- 
tem Hamiltonian Hq such that {Hj) = by definition. 

We shall assume that terms of order {H]) are small, as 
shown in Ref . , and approximate the free energy using 
only the first term of Eq. ([24l) . Neglecting the free energy 
of the bath, since it does not depend on the variational 
parameters, we find 



^B«2^c^kVk(/k-25k) 

k 

- /3"Mn[2(cosh(2/3Xc) + cosh{PEc)) 



(25) 



6 



where Ec = ■\/4Xp + A^. Minimizing Ab with respect 
to the variational parameters yields the choice 



/k = 5k 1 



A2 



sm\v{(3Ec) coth(/3wk/2) 
Wk-Bc Vcosh(2^Xc) + CQsh{PEc) 

• • (26) 

As we are interested here in the ground state (zero tem- 
perature) behaviour of the system we take the limit 
/3 — > oo to find 



/k = ffk 1 



(27) 



Having found the optimal choice for each /k in 
Eq. (j27p . we can now insert this into our expression for 
the renormalised tunneling strength, Eq. (1171) . Taking 
the continuum limit and using the same form of Ohmic 
spectral density as before, we obtain the following self- 
consistent equation 



A, 



A 



A2 



- oJcEc 



exp 



I + oJcEc 



A. (28) 



Note that with the replacement Ec f] (or 2Xc e) 
this equation is identical to Eq. ([T2|) derived in Section HIl 
when considering a single spin with finite bias. This 
stems from the observation that, from the point of view 
of one of the spins, its Ising-like coupling to the other 
spin can be thought of as providing an effective energy 
difference between its CTz eigenstates. 

The solutions of Eq. ([28]) give values of A^ that corre- 
spond to stationary points of the free energy approxima- 
tion Ab- That a given solution exists does not necessar- 
ily mean that it is appropriate to assume that the system 
will adopt this value. Rather, we assume that (within our 
approximate treatment) the system will adopt the value 
of A,, that gives the lowest Ab [44]. To see which solu- 
tion for Ar will be favored, we compute the free energy 
at zero temperature using the variational parameters we 
have just derived in Eq. 1^7^: 



A, 



-Ec - aEc 



We. 



Ec 

LUc 



(29) 



Since we are working in the limit A/lOc 1, it must also 
be true that (A^/wc)^ ^ Ec/ojc, regardless of the value 
of Xc, and we can further approximate 



Af 



Ec{l-a^ 



A2 

'e. 



(30) 



The system will adopt which ever value of A^ makes 
the second term in Eq. ([50)1 most negative, i.e. that 



value which most strongly satisfies the condition 4:X^ > 
A^(a — 1). For a < 1 it is clear that this will correspond 
to the greatest positive value of A^. Therefore, where 
multiple solutions to Eq. ([^5)) exist, for a < 1 we should 
choose the largest value of A,.. 



C. Separation-dependent localization 



In general, solving Eq. p8| for A^ analytically is not 
possible and it must be solved for numerically instead. 
However, to begin with, note that A^ = is always a 
solution, regardless of the value of a or Xc- Now, we 
can look for other analytical solutions in certain limits. 
Perhaps the simplest of these is the limit Xc — >■ 0, corre- 
sponding to two infinitely separated spins in a common 
bath or two uncoupled spins in separate baths. From 
either interpretation, we should recover the well-known 
single spin-boson results. Setting Xc = in Eq. 
gives 



exp 



-aujc 



(31) 



which in the limit A/w^, <C 1 (A^ < A) gives the well- 
established form 0, [3, [HI 



A / eA 



q/(1-q) 



(32) 



Hence, for Xc = and A/utc ^ 1, the renormalised 
tunneling strength smoothly reaches zero as a 1, and 
we predict that the critical coupling strength separating 
the delocalized and localized phases is given by ac{Xc = 
0) = 1, precisely as in the single spin-boson case. 

Let us now consider the opposite limit, Xc / A 3> 1, 
corresponding either to closely spaced spins in a common 
bath with intermediate or strong dissipation (so that Xc 
is large), or two spins in separate baths coupled via a 
relatively strong Ising interaction. Since A^ < A, we 
may also assume Xc/Ar ^ 1. Setting Ec ~ ^Xc and 
neglecting A^ in the denominators in both bracketed fac- 
tors in Eq. (^5]) , we find 



Ar 



eA2 



2u}cXc 



Q/(l-2a) 



(33) 



Within the limits this expression has been derived, the 
bracketed factor is small and we observe that A^ — 
as a — >■ 0.5. We conclude that for an Ising strength 
Xc much larger than the bare tunneling strength, the 
critical system-bath coupling strength is no longer given 
by the single spin-boson value (ac ~ 1), but instead by 
ac{Xc/A >• 1) w 0.5. This can also be seen by expand- 
ing our expression for the minimised free energy, Eq. pop , 
to lowest order in Ar/Xc, which gives 



Af 



2Xc 1 + ^(1^ 2a) 



%xl 



(34) 



> A, 



It is clear from this expression that in the limit Xc 
a finite Aj. will be favoured only if a < 0.5. 

When neither of these conditions are met, i.e. when 
Xc/A'^l, Eq. (UHl) is best studied graphically. To do 
so, we define the left hand side of Eq. (|28p as a function 
<d{Ar)- Any points at which 8(Ar) crosses the line A 
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b. 0.05 0.1 



0.00 0.02 0.04 0.06 0.08 0.10 

A. 

FIG. 3: (Color online) The renormalised tunneling strength is 
the solution to 9(Ar) = A (see Eq. ((28)) '). Main: Here we plot 
Q{Ar) for various values of a and for a fixed spin separation 
corresponding to foitd^c) = 0.05. Setting lOc = 1 then gives 
Xc = 0.025a. We see that as a is increased the solution 
for Ar- decreases until a point at which the curve B(Ar) just 
touches the line A = 0.2. After this point, the only solution 
is Ar- = 0. Inset: Here we plot the same curves but with 
foitd^c) = 0.3. In this case fo (and hence Xc) is so large 
that Qfc = 0.5, since above this value 0(Ar) never crosses A. 




FIG. 4: (Color online) Main: Numerically evaluated renor- 
malized tunneling strength as a function of the system bath 
coupling strength for various spin separations, captured by 
the value of the function /d measuring the bath correlations. 
The solid black curve corresponds to infinitely separated spins 
(no Ising interaction), fo = 0. The dashed blue curve corre- 
sponds to an intermediate spin separation (or Ising interac- 
tion), /d = 0.05. The dotted red curve corresponds to small 
spin separation (or large Ising interaction), fo ~ 0.3. Inset: 
Here we show how the critical coupling strength varies with 
the spin separation. For all curves, A/ujc ~ 0.2. 



will then give non-zero solutions for A^. In the main 
part of Fig. Owe plot Q{Ar) for a fixed spin separation, 
corresponding to a value of foitd^c) = 0.05, and for 
various values of a. Also shown is the dashed line at 0.2, 
which represents the value of the bare tunneling strength 
A (in units of ujc) taken here. The first feature to notice 
is that the curve shows a dramatic change in behavior 
as the coupling strength a moves through the value 0.5, 
developing a minimum in the first quadrant for a > 0.5. 
Therefore, when a < 0.5, we can always expect a single 
finite solution for A^. On the other hand, when a > 0.5, 
depending on the specific values of the ratio ATc/A and 
a, the curve &{Ar) may not cross the line A at all, just 
touch it, or dip low enough to cross it twice. 

In the main part of Fig. [3] the ratio Xc/A is small 
enough such that &{Ar) does dip below A for a > 0.5, 
and the critical coupling strength is then given when 
&{Ar) just touches the line A. We see that 0.5 < ac < I 
in such cases. In the inset we show the same set of 
plots, but for a smaller spin separation corresponding to 
a higher value of foitd^c) =0.3. In this case, the ratio 
Xc / A is large enough that once Q{Ar) changes its qual- 
itative behavior (i.e. when a > 0.5), it never crosses the 
line A and the only solution to Eq. is A^ = 0. This 
confirms the limiting behavior, ac ~ 0.5 for Xc/ A^ 3> 1, 
discussed earlier in reference to Eqs. (l33l) and (p4)) . 

Knowing how we expect Aj. to behave in certain limits 
we may now numerically solve Eq. (|28p to find the renor- 
malised tunneling strength as a function of a, for various 
spin separations characterized by the function fD{td^c) 
measuring the bath correlations. We shall restrict our- 
selves here to the two and three dimensional cases, and 
the results can be seen in Fig. 2] The solid black curve 



shows the renormalized tunneling for infinitely separated 
spins, i.e. /d = (no bath correlations). As expected, in 
this regime of zero Ising strength, A^ as the system- 
bath coupling strength a — 1, precisely as in the single 
spin case. The dashed blue curve shows the variation of 
the renormalized tunneling with a for a spin separation 
corresponding to Jd — 0.05. For this intermediate sep- 
aration (or, equivalently, Ising strength) we see that A^ 
discontinuously approaches zero as a reaches a critical 
value somewhere between 0.5 and 1.0 {ac ~ 0.65 for the 
values of A/wc and fo used here). This agrees with the 
intuition we gained previously from Fig. [3] The red dot- 
ted curve corresponds to a small spin separation giving 
Jd = 0.3 (or large Ising strength). Here, A^ — > contin- 
uously as a — ?> 0.5, again in agreement with our analysis 
of Fig. [3] In the inset of Fig.|3J we show explicitly how the 
critical coupling strength depends on the qubit separa- 
tion. As expected, for large spin separations [fu — 0) ac 
tends to 1 , while as the spins are brought closer together 
and fo increases, ac approaches its minimum value of 
0.5. 

The last piece of information needed to complete our 
picture is the value of /u, say /do, after which the 
crossover always occurs around Uc — 0.5 (i.e. beyond 
/do the bath correlations are large and Xc quickly dom- 
inates with increasing a). Finding where the minimum 
of 6(Ar) just crosses the line A yields the simple result 

/Do = 2e(^^^ • (^^^ 
For A/lUc = 0.2, we get /do ~ 0.22, in agreement 
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with Fig. m From this expression it can be seen that 
the further a given system hes within the scahng hmit 
(A/wc <^ 1), the larger the range of spin separations 
which differ from the single-spin case (ac ~ 1). 

D. Section summary 

We conclude this section with a brief summary. For 
distantly separated spins or negligible Ising strengths, the 
delocalized-localized crossover corresponds to the well 
studied single spin-boson model d, 0, IH, Q ■ The crit- 
ical coupling strength after which the tunneling element 
is renormalised to zero is predicted to be ac = 1, for 
A/lUc ^1. As the spins are brought closer together 
within their common bath, they become coupled via an 
Ising-like interaction. This causes the crossover region to 
drop from a,, = Ij as can be seen by tracing from left to 
right in the inset of Fig. 2] At a certain spin separation, 
the function scaling the Ising strength, foitdi-^c), reaches 
a special value, foo, given by Eq. (|35l) . For this spin 
separation, and all smaller separations, the crossover is 
predicted to occur around ac — 0.5. 

V. FULL VARIATIONAL TREATMENT 




0.0 0.2 0.4 0.6 0.8 1.0 



a 

FIG. 5: (Color online) Main: Comparison of the numerically 
calculated Ising strength (points) and the crude Ising strength 
(solid lines) as a function of a, with A/uJc = 0.2. This is done 
for tdijJc = 1 (red line, circular markers), tdCdc = 15 (blue 
line, square markers) and in the limit tduic — >■ cxd (black line, 
triangular markers, close to the a;-axis). After the system- 
bath coupling strength a — a^, Ar (see inset) and the 
crude Ising strength matches the full value. The dashed line 
shows the small Ising strength approximation of Eq. (|42p . In- 
set: Renormalized tunneling strength calculated using the full 
Ising strength (points) and crude Ising approximation (lines) 
as a function of q, for the same parameters as the main figure. 



A. Free energy minimization and self-consistent 
equations 

The results presented in the previous section were ob- 
tained by approximating the induced Ising strength, X, 
by a value Xc-, through the replacement /k — > 3k- This 
significantly simplified the task of finding the set of vari- 
ational parameters {/k}, which then allowed us to deter- 
mine the renormalised tunneling strength in a straight- 
forward manner. To go beyond this approximation, we 
shall now perform the variational calculation making no 
such simplification, and hence use the full /k-dependent 
Ising strength given in Eq. ([21]) . 

As before, we calculate the free energy associated with 
the Hamiltonian H = Hq + Hi, given by Eqs. dUl), ^ 
and (fT9|) . This leads to an expression for Ab identical to 
Eq. ((251), but with Xc replaced by X: 



Ab «2^c^kVk(/k-2.gk) 

k (36) 

- p-Hn [2(cosh(2^X) + cosh{(3E))] , 



transformed Hamiltonian H. We proceed by inserting 
Eq. (j37[) into our expressions for the renormalized tun- 
neling strength, Eq. pT| . and the full Ising strength, 
Eq. (j2ip . For simplicity, we now restrict our discussion 
to system-bath coupling in three dimensions, in which 
we may write k • (ri — r2) = ujtdC0s{9), where 6* is a po- 
lar angle in k-space over which we must integrate. We 
then obtain the following two equations which we must 
simultaneously solve self-consistcntly: 



Ar — A exp 



1 "'0 



(38) 



and 



-1 /^UJc 



X^— dx G{uj, x)(2 ~ G{uj, x)) cos{u}xtd)duj, 

(39) 

where x = cos{9), and we have defined the function 



G{uj,x) 



E + 2X cos{u}xtd) 
E + 2Xcos{ujxtd) + Alluj 



(40) 



where E — ^ AX'^ + A^. Minimization with respect 
to the variational parameters {/k} gives us the zero- 
temperature condition 

_ / i; + 2Xcos(k.(ri-r2)) \ 
^''~^\e + 2X cos(k . (n - r^)) + A^/u, J' 

which is consistent with our assumption /k = /_k, 
used with reference to Eq. (jlSp in our derivation of the 



B. Comparison of full and crude Ising strengths 



Extracting useful analytic expressions from Eqs. 
and (p9)) is not easily achieved. However, we note that the 
values Ar — and X = Xq solve these equations exactly. 
That is, in the localized regime, where the ground state 
becomes completely dominated by the Ising term, the 
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FIG. 6: (Color online) Main: Induced Ising strength as a func- 
tion of the scaled spin separation, tdiJc, for two values of the 
system-bath coupling strength, a — 0.4 (blue circular mark- 
ers) and a = 0.6 (red square markers). The markers indicate 
values calculated numerically from Eqs. (|38|) and 1)39^ and 
the solid lines represent Xc values calculated using Eq. (|22|l . 
Inset: Magnification of the lower right corner, revealing how 
the discrepancy between X and Xc increases when the Ising 
strength is small enough such that Ar 7^ 0. 



Ising strength is given by its crude value. This tallies with 
our earlier assertion that the crude Ising approximation is 
essentially a strong system-bath coupling approximation 
on the induced interaction strength. Let us also consider 
the regime in which the spins are distantly separated. 
On physical grounds, we expect that X — >• as |ri — 
— 00, since it seems inappropriate that the bath could 
mediate an interaction between spins separated by a large 
distance (certainly, we know that Xc as the spin 
separation is increased to infinity). This can be seen in 
the present case by making the assumption that for large 
td (i.e. large spin separation), X will be small (which 
we shall justify numerically in the following) and expand 
the integrand of Eq. ((39|) to second order in X. Having 
done so, the integrations with respect to x and lo can be 
performed analytically, leaving a quadratic equation for 
X which we write as 



= hoiU) + X{hi{td) - 1) + X^h2{U), 



(41) 



where /iq, /ii, and h2 are cumbersome expressions (pro- 
portional to a) which we shall not give here. Taking 
the limit td — > 00, we find that /iq — > and hi — > 
(20^^/2 (Ar -f LOcf' . Applying the same limit to /i2 is less 
straightforward, although it is easy to see graphically 
that /i2 — > as trf -> 00. Hence, as expected, we have 
confirmed that X — as — !• 00. Further, when X = 0, 
the self-consistent equation for the renormalised tunnel- 
ing strength, Eq. (|38|) . reduces to that for a single spin 
given by Eq. ([3T|l . 

When neither the spin separation nor the system-bath 
coupling strength are large enough such that the above 
arguments apply, we must solve the self-consistent equa- 
tions by numerical iteration. Solutions found in this way 
are shown in Fig. [51 where the plot points are calculated 



iteratively from Eqs. (|38)) and ([39)) . and the solid lines 
calculated using the crude Ising approximation of the pre- 
vious section. Red circular points correspond here to a 
small spin separation, i^Wc — 1, blue squares to an inter- 
mediate separation, tdUJc — 15, and black triangles to the 
limit — > cxD (for which there is no discrepancy between 
the full and crude Ising strengths). From the main part 
of the figure we can see that the crude value of the Ising 
strength generally gives a reasonably good approxima- 
tion to the full expression. As the system-bath coupling 
strength is increased, there comes a point at which the 
tunneling strength becomes entirely suppressed, — >■ 
(see figure inset), in which case Eq. p9| for X reduces 
to the simpler form of Eq. ([22]) . Hence, in the localized 
regime X = Xc, as expected. From the inset of Fig. [5] 
we see that the behaviour of the renormalized tunneling 
strength is well approximated across a range of different 
parameter regimes by replacing X by Xc in the self- 
consistent equations. Hence, our analysis of the local- 
ization crossover in the two-impurity spin-boson model 
given in the previous section is expected to hold true, 
even when the full bath-induced Ising form is used. 

In order in reproduce the behaviour of X for small val- 
ues of a and moderate spin separations, where it differs 
most markedly from Xc in Fig. [Sj we can expand the 
solution to Eq. (|4T|) to first order in a. In doing so, we 
find X w ho{td), with 



, , aAr . , . 
ho[td) = --r: — sm(tdUJc) 
2tdH 

+ :^(Ci(idM) - CiitdAr))itdAr cos(tdAr) - sm{tdAr)) 
2td 

OL 

+ — (Si(tdM) - Sl{tdAr)){tdAr Sm{tdAr) + cos{tdAr)), 

^td 



(42) 



where Ci(2;) = — cos{t)/tdt is the cosine integral 
fimction, and we have made the substitution /i — A^-t-Wc- 
The dashed curve in Fig. [5] shows this function plotted for 
tdLOc = 15, where we also approximate the renormalised 
tunneling strength as A^ ~ A(A/a;c)"^^"'^~"-'- 

Lastly, in Fig. [5] we plot a comparison of the behaviour 
oi X and Xc with varying (scaled) spin separation, tdUic- 
Recall that when A^ — X = Xc, as can be seen in 
the majority of the plot for a = 0.6. For this value of 
the system-bath coupling, a > ac over almost the full 
range of separations considered, and the tunneling is con- 
sequenctly renormalized to zero for most values of t^Wc 
too. As the spin separation is increased, X decreases, 
and there comes a point at which A^ ^ (tdUic ~ 26). 
Here, we begin to see deviations of X from Xc- When 
a = 0.4 the system is always in the delocalized regime 
(A,. 7^ 0) and we therefore see deviations of X from Xc 
for all spin separations. 



10 



VI. VARIATIONAL GROUND STATE 

A. Two-impurity Hamiltonian in the displaced 
oscillator basis 

In the preceding sections, we have used a variational 
treatment to establish how both the renormalised tun- 
neling strength and bath-induced Ising interaction vary 
as a function of system-bath coupling strength and spin- 
separation in the two-impurity spin-boson model. We 
shall now use this information to explore the interplay of 
these two quantities in determining how the form of the 
ground state of the system changes in different parameter 
regimes. From this, we shall identify a physical indicator 
of the delocalized to localized crossover in the dissipative 
two-spin system. 

To obtain the variational ground state we generalize 
the procedure given in Section to two spins. We 
write the total Hamiltonian [Eq. (fT3|) ] in a displaced 
oscillator basis, this time defined by the four states 
|00) , |01) , |10) , \B++) 111)}, with 

|S±±) =n^(±"ke''^'""')n^(±"ke^'''""')l^")' (43) 

k k 

and 

|S±t) -n^(±"ke^'''""')n^(Take^'^-""^)|i3o), (44) 

k k 

where once again = /k/'^k and \Bo) is the state of the 
bath for vanishing system-bath coupling. In this basis, 
the two-impurity spin-boson Hamiltonian becomes 

H = -^ial+al)~ 2Xdldl + 2i?, (45) 

where the zero temperature limit has been taken, and R, 
Ar, and X are defined in Eqs. ([6|), (fT7|) . and ([2T|). respec- 
tively. Diagonalizing this Hamiltonian gives a ground 
state energy of Aq = 2R — E, and corresponding ground 
state 

|<i>o) = A^oflS— )|00) + |B++)|11) 

. (46) 
-e(|i?+-)|10) + |B-+)|01))), 

where iVo = (2(1 + (^))-^/'^, ( ^ {2X - E)/Ar, and 
E = ^ylX^~{~A^ as before. Minimizing Aq with respect 
to the variational parameters leads to exactly the same 
condition [Eq. (|57|) ] as derived in Section |Vl Therefore, 
we shall make the crude Ising approximation to evaluate 
X and Ar, giving all of the required information relating 
to the variational ground state. 

B. Experimental signatures of 
localization-delocalization crossover 

To show how evidence for the localization crossover 
might be observed experimentally, in Fig. [7] we plot the 
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FIG. 7: (Color online) Expectation value of al (or a^) as a 
function of the scaled spin separation for different values of 
the system bath coupling strength a (=0.20, 0.30, 0.55, 0.65 
ordered as indicated). For a = 0.55 and a = 0.65 we see 
that at a particular spin separation there emerges a non-zero 
expectation value, signifying the crossover from localization 
to delocalization. 
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FIG. 8: (Color online) Expectation value of (or a^) as a 
function of A (measured in units of X) for a = 0.35 (black 
solid line), a = 0.6 (red dotted line), and a — 0.75 (blue 
dashed line). For these plots the induced Ising strength was 
kept at X = 0.0015ajc for each a, with ujc = 1. 



ground-state expectation value of the single-spin oper- 
ator al (or equivalently cr^), (ct^) = (<&o| ""x |*&o) ~ 
—2^{B)/{1 -I- C^), as a function of the scaled spin sep- 
aration for various values of the system-bath coupling 
strength. For small values of a (a = 0.2, 0.3) the tunnel- 
ing element is renormalised to a finite value (delocalized 
regime) and (tr^) is predominantly determined by the 
relative size of the bare tunneling element to the Ising 
strength, saturating at a value {a].) ~ {B) at large spin 
separations (small X). There is no qualitative change in 
the ground-state form as the relative size of A^. and X 
varies, in this case through increasing the spin separation. 

For larger values of a, lying between 0.5 and 1, the 
Ising strength at small spin separations is large enough 
such that the renormalized tunneling strength is com- 
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pletely suppressed, and (cr^) — (localized regime). 
As the spin separation increases, the Ising strength de- 
creases, and there conies a point at which X is small 
enough such that (cr^) can now take on non-zero values 
(delocalized) /or the same value of a. Therefore, if it 
is possible to engineer a pair of Ising-coupled spins for 
which the Ising strength can be varied, and 0.5 < a < 1, 
the crossover region should be identifiable by the emer- 
gence of a non-zero value for {a].) (or (cr^)) as the Ising 
interaction is decreased. 

It is also possible to observe the crossover behaviour 
without the need for varying the Ising strength, by in- 
stead altering the bare tunneling frequency due to the 
applied field. In Fig. |8]we again plot {a].) but this time 
as a function of the bare tunneling strength with fixed 
bath mediated Ising strength X. For a = 0.35 we ex- 
pect no crossover in ground state behaviour and we see 
[a].) — >■ only as A — >■ 0. For the curves corresponding 
to a > 0.5, when IS. /X is small we are in the regime in 
which Eq. (155)) is valid. As such, Uc — 0.5 and we see 
[a].) — 0. As the ratio A/X is increased, we eventually 
move into a regime in which Eq. ([32]) is valid and ac — > 1 . 
For a = 0.6 and a — 0.75 we must therefore enter the de- 
localized regime as A/X increases, and {a].) thus begins 
to take on non-zero values. 

C. System-bath entanglement 

Quantum phase transitions are associated with non- 
analyticity in the entanglement present in the total 
system-plus-bath state j45l - l48| . Although the variational 
treatment may not identify a true quantum phase transi- 
tion, it is expected that the change in ground-state prop- 
erties that are identified will have a manifestation in the 
entanglement [4] . Since, within the variational approach, 
the total state [Eq. (|46| ] is a pure state, we can inves- 
tigate such behaviour in our model simply by tracing 
out the bath degrees of freedom and calculating the von 
Nuemann entropy of the two-spin state. This will give a 
measure of the degree to which the spins are entangled 
with the bath [l^l ■ We define the reduced two-spin state 
as p = trB(|$o)('I'o|), where tre denotes a trace over the 
bath degrees of freedom. The von Neumann entropy is 
then defined as 

4 

5 = -pln(p) = -^Tan(r,), (47) 

i=l 

where the are the four eigenvalues of p . 

In Fig. [S] we plot the entropy as defined above (nor- 
malised by its maximum possible value) for three dif- 
ferent spin separations, corresponding to Jd = (black 
solid line), fo = 0.05 (blue dashed line) and fo — 0.3 
(red dotted line). For fo = the situation is identical 
to the single spin case. As a is increased, the extent to 
which the spins and the bath interact increases and their 
state becomes ever more entangled. For the curves cor- 
responding to fn = 0.05 and fo = 0.3 we see a similar 
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FIG. 9: (Color online) Normalized von Neumann entropy of 
the two-spin system as a function of a for /d = (black solid 
line), fo = 0.05, (blue dashed line) and fn ~ 0.3 (red dotted 
line), with A/wc = 0.2. 



situation for small values of a. However, for moderate 
values of a we see that the entanglement reaches a max- 
imum and then begins to fall. This corresponds to the 
onset of the crossover between delocalization and local- 
ization in the ground state. At the critical values of a for 
these spin separations (ac = 0.65 and Uc — 0.5, respec- 
tively, for these parameters), the entanglement sharply 
drops to a value of 0.5 as Aj. — >• 0. For a single spin 
5 = in its localized regime [J]. In the present case we 
find S'/ln4 = 0.5 since there is nothing in our model to 
lift the degeneracy between the states |00) and |11) in 
the localized regime. 



VII. SUMMARY 

We have investigated the delocalized to localized 
crossover for a pair of two-level systems in a common 
bosonic bath. Our analysis followed closely that intro- 
duced for single spins by Silbey and Harris [31] which 
used a variational approach. The crossover region is 
identified by a complete suppression of the tunneling ele- 
ment {Ar — >■ 0) as the system-bath coupling is increased 
(a ac). We find an interesting interplay between the 
magnitude of an environment-induced Ising spin interac- 
tion {X) and the applied tunneling field (A) in deter- 
mining ac- In particular, our analysis suggests that the 
presence of the Ising term encourages the spins to enter 
the localised regime at a smaller value of the system-bath 
interaction than in the single spin case. Specifically, only 
for infinitely separated spins do we recover ac = 1, as in 
the single spin-boson model. On reducing the spin sepa- 
ration from infinity, ac — >■ 0.5. Interestingly, ac reaches 
this minimum value at a finite spin separation, and re- 
tains this value for all smaller separations. We also ob- 
tained the variational ground state, and from this showed 
that a signature of the groud state crossover could be 
found in the emergence of a finite single-spin expecta- 
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tion value {(Tx) as either the spin separation or the ra- 
tio of tunnehng strength to Ising interaction is increased. 
The crossover should also be evident in the entanglement 
shared between the system and bath. 
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